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A study of charged hadron production in d+Au and Au+Au collisions is presented at various 
collision energies (^/s NJV = 19.6 to 200 GeV). Scaling and factorization features of pr and n 
distributions and vi are discussed as a function of Js NN and collision centrality. 



1 Introduction 

The study of heavy ion collisions constitutes an important part of the recent experimental 
and theoretical effort to understand the strong interaction which binds quarks and gluons into 
hadrons. Predictions based on QCD indicate the existence of a new state of matter dominated 
by the strong interaction, where bound hadrons no longer exist, provided the energy density 
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Figure 1: Ratio of pr distributions from Au+Au collisions to various reference distributions at ^/s N N =62.4 GeV 

(filled symbols) and 200 GeV (open symbols), in six centrality bins. The Raa, R A r A art , i?p^, ar * quantities are 
discussed in the text. The dashed line is the N co ii scaling expectation; the gray band is the uncertainty of N co u. 



is sufficiently high (greater than a GeV/fm 3 ). Heavy ion collisions are the only method to 
create such a high energy density in the laboratory. The current understanding of the phase 
structure of strongly interacting matter, the properties of the matter in the various phases and 
the nature of the transition between them is, to a large extent, driven by experiment. One of 
the important discoveries at the Relativistic Heavy Ion Collider (RHIC) is that an extremely 
high energy density system is created, where hadronic degrees of freedom cannot be relevant 
any more. There is evidence for a very significant level of interaction between the constituents 
of this system, as opposed to earlier expectations. 

In this experimental talk we will emphasize a few basic scaling and factorization features of 
the data collected by the PHOBOS experiment at RHIC, in comparison with earlier measure- 
ments. These simple rules highlight common features of collisions of heavy ions (Au+Au) and 
simpler systems (d+Au, p+p) in a broad range of collision energies (y/s NN =6.7 to 200 GeV). 
Details of these findings can be found in a volume summarizing the results of the four experi- 
ments from the first three years of RHIC, in which our contribution^ is titled 'The PHOBOS 
Perspective on Discoveries at RHIC. To collect the data presented here, we used the magnetic 
spectrometer and the multiplicity arrays (covering the < 5.4 pseudo-rapidity region) of the 
PHOBOS experiment (described in detail elsewhere^). 



2 Scaling features in heavy ion data 

We concentrate on three topics: transverse momentum (pr) and pseudo-rapidity (rj) distribu- 
tions of charged hadrons, and the azimuthal asymmetry of their production which is called 
elliptic flow. The consensus in the heavy ion community is that in the early stage of high energy 
heavy ion collisions a new state of strongly interacting matter has been created. The above 
observables are related to relevant physical characteristics of this newly created matter, such as 
the suppression of high-pr particles (jet quenching); the initial energy density and entropy after 
the collision and the boost-invariance of particle production along the colliding beam direction; 




Figure 2: Pseudo-rapidity density distributions of charged particles measured in p+Emulsion and d+Au collisions 
at various energies, approximately shifted to the rest frame of the lighter (a) and heavier (b) species. On the 
rightmost panel, scaling features of Au+Au data are used to predict the rj distribution at the LHC (Pb+Pb) . 

and the collective motion of particles resulting from secondary interactions and the properties 
of the equation of state. 

Determining the centrality of a heavy ion collision is essential to provide the basis of 
comparison with more elementary (p+p, p+p) processes. Instead of the impact parameter, 
two different quantities are used to quantify the centrality: the number of participant nucleon 
pairs (N p ar t/ 2) and the number of binary nucleon-nucleon collisions (N co u). Since the nucleon- 
nucleon cross section increases with collision energy, so does the N co u/N par t ratio, which also 
grows with decreasing impact parameter for simple geometrical reasons. These two quantities 
are calculated by measuring the charged hadron multiplicity in various regions of rj, combined 
with a comprehensive Monte-Carlo simulation including the Glauber model. ^ 

If one normalizes the transverse momentum distribution of charged hadrons measured 
in different centrality bins to the pt distribution measured in p+p(p) collisions and also divides 
by Ncou, one observes a gradual decrease of this quantity (Raa) with decreasing impact param- 
eter (Fig. Note that the expectation for 'hard' collisions with large momentum transfer and 
no re-interaction with the created medium would be a constant Raa equal to unity. However, 
the Raa^ quantity (where we replaced N co u by N par t in the denominator) scales with centrality 
much more precisely in Au+Au collisions at ^/s NN =62A and 200 GeV,^ pointing to new physi- 
cal interpretations.^ More strikingly, the px spectra normalized to the most central bin at each 
energy, and also normalized by N par t, as illustrated in the bottom row of Fig. ^ are identical 
within errors in each of our centrality bins. This is a clear factorization of collision energy and 
centrality. It will be interesting to compare these conclusions to the data from half a billion 
Cu+Cu events taken in the present RHIC run, since, for the same number of participants, the 
collision zone will have a very different geometry in the Cu+Cu and the Au+Au events. 

Another simple feature of the data is the extended longitudinal scaling of pseudo- 
rapidity (rj) density distributions of charged particles, in p+Emulsion and d+Au, as shown 
in Fig. [21 By plotting these distributions as a function of rf = rj ± ytarget , thus effectively looking 
at them from the rest frame of one of the colliding beams, we observe that the data at various 
energies fall on a common limiting curve, in both reference frames. The longitudinal scaling 
extends over more than an order of magnitude in beam energy. 

A similar scaling was observed earlier in heavy ion (Au+Au) collisions by PHOBOS, illus- 
trated on the right panel of Fig. |2]for the 6% most central data. In addition, applying the fact 
that dN/drj at rj ~ in central collisions scales logarithmically with ^s NN , we can extrapo- 
late the dN/drj distribution to LHC energies and give a simple experiment-based prediction at 
^5^=5500 GeV. The prediction gives dN/drj\, n= Q 1100 and 14000 charged particles in total.^ 

Longitudinal scaling has also been recently observed in the elliptic flow of particles pro- 
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Figure 3: The elliptic flow parameter V2, averaged over centrality (0-40%), as a function of 77' = r\ — ybeam (left) 
and \r)\ — ybeam (right) for various collision energies (Au+Au). Only statistical error bars are shown. 



duced in heavy ion collisions. The elliptic flow parameter (^2) is a sensitive probe of the 
properties of the newly created, very dense and hot matter in the early stage of the collision. In 
Fig. the pseudo-rapidity dependence of this V2 parameter is plotted for semi-central Au+Au 
events and for various energies, where we use the rf = rj — ybeam parameter again. There seems 
to be a universal curve governing vi over a broad range of r]'. This is shown more precisely on 
the right panel of Fig. EJ where we used the symmetry and plotted v 2 as a function of 1 77 1 — ybeam • 



3 Conclusions 



We reported a few simple scaling and factorization properties of charged hadron production 
in heavy ion collisions at RHIC energies. We observed that the transverse momentum spectra 
approximately scale with the number of participant nucleons, and the energy and centrality 
dependence of these spectra precisely factorize in the range of centrality and pt we studied. 
Longitudinal scaling was seen in an extended pseudo-rapidity region in Au+Au (and d(p)+A) 
collisions, both in case of the dN/drj pseudo-rapidity density and in case of the V2 elliptic flow 
parameter. These simple features of the data impose constraints on models attempting to 
describe and understand the basic particle production mechanism in heavy ion collisions. 
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